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£L, Abstract. We present quite a powerful method in investigations of different phenomena that 

^^ ' can appear when neutrinos and electrons propagate in background matter. This method implies 

use of exact solutions of modified Dirac equations that contain the correspondent effective 
Q\ , potentials accounting for the matter influence on particles. For several particular cases the 

exact solutions of modified Dirac and Dirac-Pauli equations for a neutrino and an electron in the 

background environment of different composition are obtained (the case of magnetized matter is 

>^— ' also considered). Neutrino reflection, trapping, neutrino pair creation and annihilation in matter 

C— . and neutrino energy quantization in a rotating medium are discussed. The neutrino Green 

' t functions in matter are also derived. The two recently proposed mechanisms of electromagnetic 

»-~H, radiation by a neutrino and an electron in matter (the spin light of neutrino and electron, SLv 

and SLe) are considered. A possibility to introduce an effective "matter-induced Lorentz force" 
acting on a neutrino and an electron is discussed. A new mechanism of electromagnetic radiation 
that can be emitted by an electron moving in the neutrino background with non-zero gradient 
of density is predicted. 
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■<^j- ■ 1. Introduction 

o. 

The problem of particles interactions under an external environment influence, provided by 
the presence of external electromagnetic fields or media, is one of the important issues of particle 
physics. In addition to possibility for better visualization of fundamental properties of particles and 
their interactions being imposed by influence of an external conditions, the interest to this problem 
is also stimulated by important applications to description of different processes in astrophysics 
and cosmology, where strong electromagnetic fields and dense matter may play an important role. 

The aim of this paper is to present a rather powerful method in investigations of different 
phenomena that can appear when neutrinos are moving in the background matter [1,2]. In 
addition, we also demonstrate how this method can be applied to electrons moving in background 
matter [3-7]. The developed new approach [4] establishes a basis for investigation of different 
phenomena which can arise when neutrinos and electrons move in dense media, including those 
peculiar for astrophysical and cosmological environments. 

The method discussed is based on the use of the modified Dirac equations for the particles wave 
functions, in which the correspondent effective potentials accounting for matter influence on the 
particles are included. It is similar to the Furry representation [8] in quantum electrodynamics, 
widely used for description of particles interactions in the presence of external electromagnetic 
fields. In this technique, the evolution operator {/fC^i,^), which determines the matrix element 
of the process, is represented in the usual form 



UF(ti,t 2 ) = Texp 



i / j^(x)A^dx 



(1) 
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where AAx) is the quantized part of the potential corresponding to the radiation field, which is 
accounted within the perturbation-series techniques. At the same time, the electron (a charged 
particle) current is represented in the form 

i4z) = f[*e7/,,*e], (2) 

where ^e are the exact solutions of the Dirac equation for the electron in the presence of external 
electromagnetic field given by the classical non-quantized potential A^ (x) : 

{y(tfy-e^f(aO)-m e }tte(a;) = 0. (3) 

Note that within this approach the interaction of charged particles with the external 
electromagnetic field is taken into account exactly while the radiation field is allowed for by 
perturbation-series expansion techniques. A detailed discussion of the use of this method can be 
found in [9] . Many processes with electrons under the influence of external electromagnetic fields 
were investigated using this method. In particular, this method was applied [10] for derivation 
of an electron dispersion relation in external electromagnetic fields as well as in studies of the 
problem of the electron anomalous magnetic moment in external fields (see [11] for a review). 

In Section 2.1 we derive the modified Dirac equation for the neutrino wave function in the 
presence of matter and find its exact solutions including the neutrino energy spectrum (Section 
2.2). On this basis we discuss the neutrino reflection, trapping and also neutrino pair annihilation 
and creation in matter (Section 2.3). In Section 2.4 we consider the modified Dirac equation for the 
case when neutrino propagates in rotating matter and find that its energy is quantized very much 
similar to the electron energy Landau quantization in a magnetic field. In Section 2.5 we consider 
the Dirac-Pauli equation for a neutrino moving in matter. The correspondent neutrino energy 
spectrum, as well as the one of the modified Dirac equation, can be used for obtaining the correct 
values for the flavour and hclicity neutrino energy difference in matter (Section 2.6). In Section 
2.7 we use the modified Dirac-Pauli equation to get neutrino energy spectrum in magnetized 
and polarized matter. Section 2.8 is devoted to discussion of the modified Dirac equation for 
a Majorana neutrino in matter. Neutrino Green functions, for both Dirac and Majorana cases, 
are derived in Section 3. In Section 4 we apply the developed method of exact solutions of the 
quantum wave equation to the study of an electron moving in background matter and found exact 
solutions of the correspondent Dirac equation. In Section 5 we illustrate how the obtained exact 
solutions can be used in studies of different processes in matter. As two examples, we discuss 
evaluation of quantum theory of the spin light of neutrino (SLi>) and spin light of electron (SLe) 
in matter, the two recently discussed new mechanisms of electromagnetic radiation produced by a 
neutrino and an electron moving in matter. A possibility to introduce an effective "matter induced 
Lorentz force" acting on a neutrino and an electron is discussed in conclusions Section 6. We also 
predict a new mechanism of electromagnetic radiation that can be emitted by an electron moving 
in the neutrino background with non-zero gradient of density. The proposed mechanism of the 
electromagnetic radiation can be important in physics of neutron stars, gamma-ray bursts and 
black holes. 

2. Quantum equations for neutrino in matter 

2.1. Modified Dirac equation for neutrino in matter 

In [1] (see also [2,3]) we derived the modified Dirac equation for neutrino wave function exactly 
accounting for the neutrino interaction with matter. Let us consider the case of matter composed of 
electrons, neutrons, and protons and also suppose that the neutrino interaction with background 
particles is given by the standard model supplied with the singlet right-handed neutrino. The 
corresponding addition to the neutrino effective interaction Lagrangian is given by 

AL e// = -/"(*7 M ±±2!i/), ,r = V2G F ]T tfqf+^qf, (4) 

f=e,p,n 
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where 

9? = (4? - 2QO sin 2 ^ + ft,), g f = -(1$ + *„), <5 e/ = { J £ J^ p (5) 

Here I^ L and Q^ are, respectively, values of the isospin third components and the electric charges 
of matter particles (/ = e,n,p). The corresponding currents j* and polarization vectors A^ are 

% = (n f) n f vf), A? = f n f (CfVf),nfCfJl^ f + n t y t^A , (6) 

V 1 + ^Jl-vjJ 

where 6w is the Weinberg angle. In the above formulas ([6]), n/, V/ and £/ (0 < |C/| 2 < 1) stand, 
respectively, for the invariant number densities, average speeds and polarization vectors of the 
matter components. Using the standard model Lagrangian with the extra term ((4]), we derive the 
modified Dirac equation for the neutrino wave function in matter [1]: 

{ V - \l,(l + 7b)/" - m}*(x) = 0. (7) 

This is the most general form of the equation for the neutrino wave function in which the 
effective potential V^ — \{1 + 75)/^ includes both the neutral and charged current interactions of 
neutrino with the background particles and which can also account for effects of matter motion 
and polarization. It should be mentioned that other modifications of the Dirac equation were 
previously used in [12-18] for studies of the neutrino dispersion relations, neutrino mass generation 
and neutrino oscillations in the presence of matter. 

2.2. Neutrino quantum states in matter 

In the further discussion below we consider the case when matter is compose of electrons 
and no electromagnetic field is present in the background. We also suppose that the matter is 
unpolarizcd, A M = 0. Therefore, the term describing the neutrino interaction with the matter is 
given by 

/" = £|(n,nv), (8) 

where we use the notation Gf = G\f(1 + 4 sin 2 0w)- 

For the stationary states of the equation ([7]) we get [1] 

y(r,t) = e-^ t -rtu(p,E 6 ), (9) 

where u(p, E £ ) is independent on the coordinates and time. Upon the condition that the equation 
([7]) has a non-trivial solution, we arrive to the energy spectrum of a neutrino moving in the 
background matter: 



,21- ™ X2 



E e = £m/p 2 ( 1 — sa — ) + m 2 + am, (10) 

where we use the notation 

and also introduce the value n =sign(l — sa— ) in order to provide a proper behavior of the wave 
function in the hypothetical massless case. The values s = ±1 specify the two neutrino helicity 
states, v + and z/_. In the relativistic limit the negative-helicity neutrino state is dominated by the 
left-handed chiral state (z/_ « v£), whereas the positive-helicity state is dominated by the right- 
handed chiral state (1/+ rs z^.The quantity e — ±1 splits the solutions into the two branches that 
in the limit of the vanishing matter density, a — » 0, reproduce the positive- and negative- frequency 
solutions, respectively. It is also important to note that the neutrino energy in the background 
matter depends on the state of the neutrino longitudinal polarization, i.e. in the relativistic case 
the left-handed and right-handed neutrinos with equal momenta have different energies. 
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We get the exact solution of the modified Dirac equation in the form [1] 

yl + E e m am y 



p -i(E c t-pr) 

* E ,P, s (r,t) = 



2Li 



„P3 
S P 



i - ^— , 1 - sf e u 



serf 



,iii 



(12) 



V "V 1 - ^^ v /rT ^ e " / 



where the energy 2? e is given by (|10p . L is the normalization length, and 6 — arctanp2/pi- In 
the limit of vanishing density of matter, when a — ► 0, the wave function (|12j) transforms to the 
vacuum solution of the Dirac equation. 

Let us now consider in some detail properties of a neutrino energy spectrum (|10p in the 
background matter that are very important for understanding of the mechanism of the neutrino 
spin light phenomena. For the fixed magnitude of the neutrino momentum p there are two values 
for the "positive sign" (e = +1) energies 



E s=+1 = . p 2 1 -a— ) + m 2 + um, E 3 ^ 1 = Jp 2 ( 1 + a— ) + m 2 + am, (13) 





that determine the positive- and negative-helicity eigenstates, respectively. The energies in (fH 
correspond to the particle (neutrino) solutions in the background matter. The two other values 
for the energy, corresponding to the negative sign e = — 1, are for the antiparticle solutions. As 
usual, by changing the sign of energy, we obtain the values 



E s=+1 = Jp 2 (l-a—\ +m 2 -am, E s =- 1 = Jp 2 (l + a—J + m 2 - am, (14) 

that correspond to the positive- and negative-helicity antineutrino states in the matter. The 
neutrino dispersion relations in matter exhibits a very fascinating feature (see also [15, 16]): the 
neutrino energy may has a minimum at non-zero momentum. It may also happen that the 
neutrino group and phase velocities are oppositely directed. The expressions in (|13p and (|14p 
would reproduce the neutrino dispersion relations of [16], if the contribution of the neutral-current 
interaction to the neutrino potential were omitted. 

In the general case of matter composed of electrons, neutrons and protons the matter density 
parameter a for different neutrino species is 

ttv e ,v u ,u r = — 7= (n e (4sin 2 6W + O) +n p (l - 4sin 2 9 W ) - n n ) , (15) 

2V2 m V ) 

where g = 1 for the electron neutrino and g = — 1 for the muon and tau neutrinos. 

Note that on the basis of the obtained energy spectrum (|10p the neutrino trapping and 
reflection, the neutrino-antineutrino pair annihilation and creation in a medium can be studied 
[15,16,19-21]. 

2.3. Neutrino reflection, trapping and neutrino-antineutrino pair annihilation and creation in 
matter 

Analysis of the obtained energy spectrum (fT"3|) . (fT4| enables us to predict some interesting 
phenomena that may appear at the interface of the two media with different densities and, in 
particular, at the interface between matter and vacuum. Indeed, as it follows from (|13p and (|14p 
(see also [20]), the band-gap for neutrino and antineutrino in matter is displaced with respect to 
the vacuum case in neutrino mass and is determined by the condition am — m < E < am + m. 
For instance, if a — «2 > 2 then there is no band-gap overlapping. This situation is illustrated in 
Fig.l. 

Let us consider first a neutrino moving in the vacuum towards the interface with energy 
that falls into the band-gap region in matter. In this case the neutrino has no chance to survive 
in the matter and thus it is reflected from the interface. The same situation is realized for the 
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+m 
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Figure 1. The interface between the vacuum (left-hand side of the picture) and the 
matter (right-hand side of the picture) with the corresponding neutrino band-gaps shown. 
The parameter a = a^ > 2. 



antincutrino moving in the matter with energy falling into the band-gap in the vacuum. In this 
case the antineutrino is trapped by the matter. When the energies of neutrino in the vacuum 
or antineutrino in the medium fall into the region between the two band-gaps the effects of the 
neutrino-antincutrino annihilation or pair creation may occur (see the first paper of Rcf. [16] 
and [19-21]). Indeed, the "negative sign" energy levels in the matter (the right-hand side of 
Fig.l) have their counterparts in the "positive sign" energy levels in the vacuum (the left-hand 
side of Fig.l). The neutrino-antineutrino pair creation can be interpreted as a process a particle 
state appearance in the "positive sign" energy range accompanied by appearance of the hole state 
in the "negative sign" energy sea. The phenomenon of neutrino-antineutrino pair creation in the 
presence of matter is similar to the spontaneous electron-positron pair creation the electrodynamics 
according (Klein's paradox). 



2.4- Neutrino quantum states in rotating medium 

In this section we apply our method to a particular case when neutrino is propagating in a 
rotating medium of constant density [22] . Suppose that a neutrino is propagating perpendicular to 
uniformly rotating matter composed of neutrons. This can be considered for modelling of neutrino 
propagation inside a rotating neutron star. The corresponding modified Dirac equation for the 
neutrino wave function is given by ([7]) with the matter potential accounting for rotation, 

r = -G(n,nv), v=(w»,0,0), (16) 

Here u> is the angular frequency of matter rotation around OZ axis, it also is 



where G 



G 

V2 



(17) 



supposed that the neutrino propagates along OY axis. For the neutrino wave function components 
$(x) we get the from the modified Dirac equation ([7]) a set of equation^, 

[i (d - d 3 ) + Gn]$i + [- (idi + d 2 ) + Gnu)y]^ 2 = m^ 3 , 

[ (-ifli + d 2 ) + Gnuiy] *i + [i (d + d 3 ) + Gn] W 2 = m# 4 , 

i(do + d 3 )^ 3 + (idi + d2)^4 = m^i, 

(idi - d 2 ) *3 + i (do - d 3 ) *4 = m# 2 - 

In general case, it is not a trivial task to find solutions of this set of equations. 

The problem is reasonably simplified in the limit of a very small neutrino mass, i.e. when the 
neutrino mass can be ignored in the left-hand side of (JX7J) in respect to the kinetic and interaction 
terms in the right-hand sides of these equations. In this case two pairs of the neutrino wave 
function components decouple one from each other and four equations (jTTJ) disintegrate to the two 
independent sets of two equations, that couple together the neutrino wave function components 
in pairs, (*i, W 2 ) and (* 3 , * 4 ). 



| The chiral representation for Dirac matrixes is used. 
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The second pair of equations (|17[) does not contain a matter term and is attributed to the 
sterile right-handed chiral neutrino state, ^^. The corresponding solution can be taken in the 
plain-wave form 

$ R ~ L~2 exp{i(-p t + pxx + p 2 y + p^z)}^, (18) 

where p^ is the neutrino momentum. Then for the components ^3 and ^4 we obtain from (|17p 
the following equations 

(PO - P3) *3 - (Pi - IP2) *4 = 0, 
- (pi + ip 2 ) *3 + Oo + Ps) *4 = 0. 

Finally, from (fl9|) for the sterile right-handed neutrino we get 

/ ° \ 



(19) 



#/ 



Ly^2p ( Po - P3 ) 



(20) 





Pi + 1P2 

V P3-P0 J 

(t,x,y, z). This solution, as it should be, has the 



where px = p^ 1 , p^ = (po,Pi,P2,P3) and x M 
vacuum dispersion relation. 

In the neutrino mass vanishing limit the first pair of equations (| 1T|) corresponds to the active 
left-handed neutrino. The form of these equations is similar to the correspondent equations for a 
charged particle (e.g., an electron) moving in a constant magnetic field B given by the potential 
A = (By, 0, 0) (see, for instance, [9]). To display the analogy, we note that in our case the matter 
current components nv plays the role of the vector potential A. The existed analogy between 
an electron dynamics in an external electromagnetic field and a neutrino dynamics in background 
matter is further discussed in the conclusion (Section 6). 

The solution of the first pair of equations (Tl7|) can be taken in the form 

* L - - exp{i(-p t + pix + P3z)}ip(y), (21) 

and for the components ^i and ^ 2 of the neutrino wave function we obtain from (|17[) the following 
equations 



(22) 



Po+P3 + Gnp 1 - y /p[^ j -r ] )^2 = 0, 



Vp{£-, 



V *! + (po - p 3 + Gn W 2 = 0, 



where 



V = VP [X2 + 



El 
P 



p = Gnio. 



For the wave function we finally get 



*; 



/9-*e 



-ipot+ip 1 x+ip 3 z 



V(Po -P3 + Gn) 2 + 2pN 



( (po- p3 + Gn)u N (n) ^ 

-^/2p~Nu N - X (ri) 







(23) 



(24) 



where un{t]) are Hermite functions of order N . For the energy of the active left-handed neutrino 
we get 



Po 



Pi 



2pN-Gn, N = 0,1,2,.. 



(25) 



The energy depends on the neutrino momentum component p 3 along the rotation axis of matter 
and the quantum number N that determines the magnitude of the neutrino momentum in the 
orthogonal plane. For description of antineutrinos one has to consider the "negative sign" energy 
eigeinvalues (see similar discussion in Section 2.2). Thus, the energy of an electron antineutrino 
in the rotating matter composed of neutrons is given by 



Po 



Pi 



2pN + Gn, N = 0,1,2,... . 



(26) 
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Obviously, generalization for different other neutrino flavours and matter composition is just 
straightforward (see (O and (US])). 

Thus, it is shown [22] that the transversal motion of an active neutrino and antineutrino is 
quantized in moving matter very much like an electron energy is quantized in a constant magnetic 
field that corresponds to the relativistic form of the Landau energy levels (see, for instance, the 
first book of [9]). Consider again antineutrino. The transversal motion momentum of is given by 



p± = y/2pN. (27) 

The quantum number N determines also the radius of the antineutrino quasi-classical orbit in 
matter (it is supposed that N 3> 1 and p% = 0) , 



It follows that antineutrinos can have bound orbits inside a rotating star. To make an estimation 
of magnitudes, let us consider a model of a rotating neutron star with radius Rns — 10 km, matter 
density n = 10 37 cm~ 3 and angular frequency u> — 2n x 10 3 s _1 . For this set of parameters, the 
radius of an antineutrino orbits is less than the typical star radius Rns if the quantum number 
N < N max = 10 10 . Therefore, antineutrinos that occupy orbits with N < 10 10 can be bounded 
inside the star. The scale of the bounded antineutrinos energy estimated by (f2"6")l is of the order 
p~o ~ 1 eV. It should be underlined that within the quasi-classical approach the neutrino binding 
on circular orbits is due to an effective force that is orthogonal to the particle speed. Note that 
there is another mechanism of neutrinos binding inside a neutron star when the effect is produced 
by a gradient of the matter density [19] (see also the conclusion). 

2. 5. Modified Dirac-Pauli equation for neutrino in matter 

To derive the quantum equation for a neutrino wave function in the background matter we 
start with the well-known Dirac-Pauli equation for a neutral fermion with non-zero magnetic 
moment. For a massive neutrino moving in an electromagnetic field F^ v this equation is given by 

V 

where m and /i are the neutrino mass and magnetic moment [23] [§|, a iiv — i/2(j fJ/ j u — 7^7^). 
It worth to be noted here that Eq. (|2"9]) can be obtained in the linear approximation over the 
electromagnetic field from the Dirac-Schwinger equation, which in the case of the neutrino takes 
the following form [25]: 

{ij^df, - m)*(x) = / M F (x', x)$(x')dx', (30) 



17"^ - m - -a^F^ )*(») = 0, (29) 



where M F (x',x) is the neutrino mass operator in the presence of the external electromagnetic 
field. 

Recently in a series of our papers [26,27] (see also [28]) we have developed the quasi-classical 
approach to the massive neutrino spin evolution in the presence of external fields and background 
matter. In particular, we have shown that the well known Bargmann-Michel-Telegdi (BMT) 
equation [29] of the electrodynamics can be generalized for the case of a neutrino moving in the 
background matter and being also under the influence of external electromagnetic fields. The 
proposed new equation for a neutrino, which simultaneously accounts for the electromagnetic 
interaction with external fields and also for the weak interaction with particles of the background 
matter, was obtained from the BMT equation by the following substitution of the electromagnetic 
field tensor F M „ = (E,B): 

§ For the recent studies of a massive neutrino electromagnetic properties, including discussion on the neutrino 
magnetic moment, see Ref. [24] 
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where the tensor G^ u = (— P,M) accounts for the neutrino interactions with particles of the 
environment. The substitution (|3ip implies that in the presence of matter the magnetic B and 
electric E fields are shifted by the vectors M and P, respectively: 

B^B + M, E^E-P. (32) 

We have also shown [26, 27] how to construct the tensor Guv with the use of the neutrino speed, 
matter speed, and matter polarization 4-vectors. 

Now let us consider the case of a neutrino moving in matter without any electromagnetic field 
in the background. Starting from the Dirac-Pauli equation (|29[) for a neutrino in electromagnetic 
field F UVl we apply the substitution (|31[) which now becomes 

As a result of this substitution, we obtain the quantum equation for the neutrino wave function 
in the presence of the background matter in the form [30] 



(*7"0 P 



m-^G^ )*(.,-)=(). 



(34) 



that can be regarded as the modified Dirac-Pauli equation. 

Consider an explicit solution of the obtained equation (|34|) for the case of an unpolarized 
matter composed of only electrons we have 

( \ 



G nu = 



Gi 



2y/2fi 



"fii 



V 






-ft 02 





03 -0 X 





-02 01 



, 7 = (i-/3 2 r 1/2 , 



(35) 



/ 



where (3 = (0i,02, 03) is the neutrino three-dimensional speed and n denotes the number density 
of the background electrons. From (|3"B"|) and two equations, (|29|) and (fM)) . it is possible to see 
that the term Z jn/3 in (|34|) plays the role of the magnetic field B in ((29|) . The corresponding 
neutrino energy spectrum is 



E = v p 2 (l + on 2 ) + m? — 2amps, 
This expression can be transformed to the form 



a 



— i=Gf — • 
2V2 m 



(36) 



E = 



ap\- 



'p z + m'\L - s-^-j , (37) 

that can be obtained from the neutrino vacuum spectrum by the formal shift of the neutrino mass 
m — ► ml 1 — s^2 ) . 

can be obtained in the following form [30] : 



The exact solution of the Dirac-Pauli equation 

/ 



*P,.(',*) 



i(Bt-pr) 

~2l1 



1 



1 



m— sap 
E 



m—sap 
E 






l + *a 



V 



-| in— sap 

1 S 



1 



S P3 e M 



(38) 



/ 



0, this wave function transforms to the vacuum 



In the limit of vanishing matter density, when a 
solution of the Dirac equation. 

The obtained neutrino energy spectrum (1361) . for not extremely high matter densities a^Si ^ 

1, yields the correct result for the energy difference AE = E(s = —1) — E(s = +1) of the two 
neutrino helicity states: 



AE 



o P 

E 



(39) 



where we use the notation E$ = yp 2 + m 2 . Therefore, on the basis of the obtained exact solution 
for the neutrino wave function in the case of relativistic neutrinos one can derive the probability 
of spin oscillations vl ^-> vr in transversal magnetic field with the correct form of the matter 
term [31]. 
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2.6. Flavour and helicity neutrino energy difference in matter 

Although the neutrino energy spectra corespondent to the modified Dirac and Dirac-Pauli 
equations, (|7|) and (f3~4"|) , are not the same, an equal result given by (|5§|) for the energy difference 
AE = E(s = — 1) — E(s = +1) of the two neutrino helicity states can be obtained from both of 

the spectra in the low matter density or high energy limit a^r <C 1. 


It should be also noted that for the relativistic neutrinos the energy spectrum for the neutrino 

helicity states of Eq. (|10p in the low density limit reproduces the correct energy values for the 

neutrino left-handed and right-handed chiral states: 

E VL a E(s = -1) « E a + -j=n, E VR « E(s = -1) » E , (40) 

as it should be for the active left-handed and sterile right-handed neutrino in matter. 

We should like to note, that the obtained spectra for the flavour neutrinos of different helicities 
in the presence of matter enables one to reproduce the well-known result for the energy difference 
of two flavour neutrinos in matter. In order to demonstrate this we expand the expressions for 
the relativistic electron and muon neutrino energies (given by (|10p for the Dirac case or by (|53[) 
below for the Majorana case), over m/p <C 1 and get 

E°=-^~E Q + 2a u ^m. (41) 

Then the energy difference for the two active flavour neutrinos is 

AE = E s =~ x - E a =- X = V2G F n e . (42) 

Analogously, considering the spin-flavour oscillations v &L <^ v ^ R i for the corresponding energy 
difference we find: 

AE = ET c - x - El= +l = V2G F (n e - ^n n ). (43) 

These equations enable one to get the expressions for the neutrino flavour and spin-flavour 
oscillation probabilities with resonance dependence on matter density in the complete agreement 
with the results of [31,32]. 

2. 7. Modified Dirac-Pauli equation in magnetized and polarized matter 

It is also possible to generalize the Dirac-Pauli equation (|29|) (or (j34|) ) for the case when 
a neutrino is moving in a magnetized background matter. For this case (i.e., when the effects 
of matter and magnetic field on neutrino have to be accounted for simultaneously) the modified 
Dirac-Pauli equation is [30] 

{iVfl/i -m- ^(F„u + Gp,)}*(x) = 0. (44) 

The neutrino energy in the magnetized matter can be obtained from (|36|) by the following 
redefinition 

a^a' = a-\ % (45) 

P 

where B\\ = (Bp)/p is the longitudinal to the neutrino momentum magnetic field component. 
Thus, the neutrino energy in this case reads 



E= Jp 2 + m 2 (l-s aP + fiBli y. (46) 

For the relativistic neutrinos the expression of Eq. (|46|) gives, in the linear approximation over 
the matter density and the magnetic field strength, the correct value (see [26-28]) for the energy 
difference of the two opposite helicity states in the magnetized matter: 

A e// = ^|n + 2^!. (47) 
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Note that the problem of the neutrino dispersion relation in an external magnetic field and 
matter was also studied previously in many papers with use of different methods [34] . 

Now we can consider the neutrino spin oscillations in the presence of non-moving matter being 
under the influence of an arbitrary constant magnetic field B = B +Bj_, here Bj^ is the transversal 
to the neutrino momentum component of the external field. In the adiabatic approximation the 
probability of the oscillations vl *-* vr can be written in the form, 



Pvl-*vr ( x ) = sm2 2 ^e// sin 2 



J eff 



20 eff ~ 



E hf 



K 2 



A 2 ' 



L eff ~ 



2vr 



E eff 



A 2 



(48) 



2fiB± (terms 



7 



are omitted here), and x is the distance traveled by the 



where E e ff 
neutrino. 

Let us now shortly discuss the effect of matter polarization. Consider the case of matter 
composed of electrons in the presence of such strong background magnetic field so that the 

2 

following condition is valid B > ^|-, where pf = \//i 2 — to 2 , [i and m e are, respectively, the 
Fermi momentum, chemical potential, and mass of electrons. Then all of the electrons occupy the 
lowest Landau level, therefore the matter is completely polarized in the direction opposite to the 
unit vector ||. From the general expression for the tensor G M „ (see the second paper of [26]) we 
get [30] 
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(49) 



Thus, the modified Dirac-Pauli equation (|4"4"]) with the tensor G^ given by (|4"9"]) can be used for 
description of the neutrino motion in matter which is magnetized and totally polarized in respect 
to the magnetic field vector B direction. The neutrino energy in such a case can be obtained from 
(f3"6"| by the following redefinition 



• I B \\ 
sign {-£ 

1 + sin 2 40 



ir 



lxB\\ 
P 



(50) 



In Eq. ([50|) . the second term in brackets accounts for the effect of the matter polarization. It follows, 
that the effect of the matter polarization can reasonably change the total matter contribution to 
the neutrino energy (|4l)|) (see also [33]). 

Note that the problem of the neutrino dispersion relation in an external magnetic field and 
matter was previously also studied in many papers with use of different methods [34] . 



2.8. Majorana neutrino 

We have considered so far the case of the Dirac neutrino. Now let us turn to the Majorana 
neutrino [20]. For a Majorana neutrino we derive the following contribution to the effective 
Lagrangian accounting for the interaction with the background medium 

AL eff = -riv-frfv), (51) 

which leads to the Dirac equation 



{V - 7//75.P - m}*(x) 



0. 



(52) 



This equation differs from the one, obtained in the Dirac case, by doubling of the interaction term 
and lack of the vector part. The corresponding energy spectrum for the equation (|52p is: 



E F 



y i 



2sa— 
P 



(53) 



From this expression it is clear, that the energy of the Majorana neutrino has its minimal value 
equal to the neutrino mass, E = to. This means that no effects are anticipated for the Majorana 
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neutrino such as the Dirac neutrino has at the two media interface and which are discussed above. 
So that, in particular, there is no Majorana neutrino trapping and reflection by matter. It should 
be noted that the equation (|52p and the Majorana neutrino spectrum in matter were discussed 
previously also in [16,35]. 

3. Neutrino Green function in matter 

The neutrino Green function, along with the wave function, is an important characteristic of 
the neutrino (propagation) in matter. Developing further the method of the exact solutions for 
the studies of the neutrino propagation in matter, we consider explicit Green functions for the the 
modified Dirac equation for the Dirac and Majorana neutrinos [36]. For the Dirac and Majorana 
neutrino Green functions we obtain the same equations as for the correspondent wave functions, 
([7]) and (15"2")) , with the only difference that —5(x) functions stay on the right hand sides. In the 
momentum representation the equation for the Green function has the following form: 

{iJ^-^(a + l5 )F-m}G(p) = ~l , (54) 

where a = 1 for the Dirac neutrino and a = for the Majorana case. Squaring the left hand side 
operator, it is possible to obtain the following expression for the Green function of neutrino in 

matter: 

- (q 2 - TO 2 ) (q + m) + f(q- to) P l (q + to) - f 2 qP L + 2 (fq) P R (q + to) 

&(q) = — s ! (55) 

(q 2 — to 2 ) — 2 (fq) (q 2 — to 2 ) + f 2 q 2 

where 

q=p-\(a- 1)/, q = q M >f, ? = <l»<f \ / = /m7 M , (/«) - U<f, Pl,R = \(l ± 7s) • (56) 

Now let us consider the denominator of expression (|55p. The poles of the Green function 
determine the neutrino dispersion relation. Equating the denominator to zero, we obtain quadratic 
equation relative to q$: 

(q 2 -m 2 ) 2 -2(fq)(q 2 -m 2 )+f 2 q 2 = Q. (57) 

In some special cases equation (J5TJ) can be solved analytically. One of such cases is that of uniform 
medium, moving at constant speed v parallel to the neutrino momentum p. In this case we can 
solve equation (|57|) for go and then find pa, 



1 

Po=2 



o/o + 8 | f | + e J Am 2 + (2| p - \(a - l)f | - s(f + s\f\ 



(58) 



There are four solutions of ([55]) corespondent to s = ±1 and e = ±1. From Eg.([55|) one can find, 
that all solutions except one are of definite sign for any |p|. The sign of po for e = — 1 and s = 1 
however can be both positive and negative for different |p|. One can also note, that in case of 

o/o + | f | < 2m, (59) 

the sign of this po is always negative. 

In case condition Eq. (|5"9")) holds, the solution of equations ([7]) and (|B"2"|) can be expressed in 
the form of the superposition of plane waves each with definite sign of energy. Note that if the 
condition (|59|) is violated then there exists a plane wave that has positive energy for some |p| 
and negative for others. Stated in other words, the condition (|59[) means that Green function 
(|55p can be chosen causal by imposing special rules of poles bypassing (negative poles should be 
bypassed from below and positive poles should be bypassed from above). Once we got the causal 
Green function the perturbation technique can be developed for the description of the neutrino 
propagation in matter. 

Another way to interpret the condition (|59[) is to turn attention to [20] where it was shown, 
that for the matter at rest, the spontaneous vv pair creation can take place only when /o > 2m. 
From the analysis of the allowed energy zones for neutrino in matter it follows that vv pair creation 



Method of wave equations exact solutions in matter.. 



12 



in moving matter can take place only when a/o + |f | > 2m. So that the possibility of using the 
neutrino Green function (|55[) is limited by the particular value of matter density when vv pair 
creation processes become available. 

For the Majorana neutrino moving through uniform matter at rest the condition (|59p is always 
valid for any matter densities /o because a = and f = in this case. 

4. Electron wave function and energy spectrum in matter 



In [3,5-7], it has been shown how the approach, developed at first for description of a neutrino 
motion in the background matter, can be spread for the case of an electron propagating in matter. 

Let us consider an electron having the standard model interactions with particles of electrically 
neutral matter composed of neutrons, electrons and protons. This can be used for modelling a 
real situation existed, for instance, when electrons move in different astrophysical environments. 
We suppose that there is a macroscopic amount of the background particles in the scale of an 
electron de Broglie wave length. To further simplify the model, we consider the case of nuclear 
matter [21,37] composed of neutrons. Then the addition to the electron effective interaction 
Lagrangian is 



alw = -r 



eitJ.- 



l-4sin 2 <V + 7 5 



(60) 



where the explicit form of / M depends on the background particles density, speed and polarization 
and is determined by ([4J and ([5]). The modified Dirac equation for the electron wave function in 
matter is [3] 



{*7/*0" - 2^(! - 4sin 2 6 W + l5 )f» - m e }* e (x) 



= 0, 



where for the case of an electron moving in the background of neutrons 

Gp 



f" 



-/'' 



V2 



(?' M '- A M ). 



We consider below unpolarized neutrons so that 
Gp 



f" 



v/2 



(n„,n„v) 



(61) 



(62) 



(63) 



here n n is the neutrons number density and v is the speed of the reference frame in which the 
mean momentum of the neutrons is zero. The corresponding electron energy spectrum in the case 
of unpolarized matter at rest is given by 



Ei e *> = 




+ ca n m e , a n 



1 



Lrp , 



2V% m 

where c = 1 — 4 sin 9w> Ve =sign(l — sa„— ), m e and p are the electron mass anc 
For the wave function of the electron moving in nuclear matter we get [5-7] 

/ 



(64) 
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(65) 



/ 



The exact solutions of this equation open a new method for investigation of different quantum 
processes which can appear when electrons propagate in matter. 
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5. Neutrino and electron spin light in matter 

In this section we illustrate how the developed method based on the use of the exact solutions 
of the modified Dirac equations for neutrino and electron wave functions can be used in the study 
of different phenomena that arris when a neutrino or electron move in matter. As an example, 
we discuss below the spin light of neutrino (SLv) and the spin light of electron (SLe), new types 
of electromagnetic radiation that can be produced by the Dirac particles while moving in the 
background matter. The spin light of neutrino in matter, one of the four new phenomena studied 
in our recent papers (see for a short review [38]), is an electromagnetic radiation that can be 
emitted by a massive neutrino (due to its non-zero magnetic moment) when the particle moves 
in the background matter. Within the quasi-classical treatment the existence of the SLv was 
first proposed and studied in [39] , while the quantum theory of this phenomenon was developed 
in [1-3,20,28,30,40]. The spin light of electron in matter [3,5-7] also originates from the particle 
magnetic moment procession in matter. Note that the term "spin light of electron" was used first 
in [41] for designation of the particular spin-dependent contribution to the electron synchrotron 
radiation power. It should be stressed that the SLv and SLe in matter are really new mechanisms 
of electromagnetic radiation of quite a different nature that ones considered before including the 
Cherenkov radiation of particles in medium. In particular, the spin light processes may proceed 
even when the photon refractive index in matter equals to n 7 = 1 . 

The corresponding Feynman diagram of these processes is shown in Fig. 2. The particles initial 




Figure 2. The SLv and SLe radiation diagram. 

ipi and final ipf states (shown by "broad lines") are exact solutions of the corresponding Dirac 
equations for the neutrino and electron in matter that account for the particles interaction with 
matter. The amplitude of the SLv process is given by 

/ p ihx 
J 4 ™j, f M(t P *) ^(x), f = iuj{ [S x a] + i7 5 S}, (66) 

V2luL 3 

where \x is the neutrino magnetic moment, fc M = (to, k) and e* are the photon momentum and 
polarization vectors, it — k/w is the unit vector pointing in the possible direction of the emitted 
photon propagation. The amplitude of the process SLe is given by 



I k:r 



S }1 = -W47T / dM/(z)(7^*)^==V>i(z), (67) 

where — e is the electron charge. The further evaluation of the SLv and SLe characteristics of 
the processes, such as the differential and total rates and powers, angular distributions etc, can 
be found in the mentioned above papers. From the energy momentum conservation in the SLv 
and SLe processes we obtain the following values for the spin light radiation energy 

2a n mip [E - (p + a n mi) cos Qsl] g, „ 1 „ n n 

usl = — -h — — ~ — : ~ — , E = E-ca n m h a n = — j=G F — , (68 

(E -p cos 6sl) 2 - {a n mi) 2 2^ mi 

where 6sl is the angle between possible directions of the radiation and the initial particle 
momentum p, for the case of neutrinos mi = m v and cj = c„ = 1, whereas for electrons mi = m e 
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and cj = c e = 1 — 4 sin 2 0^ . From ([55)1 it follows that for the relativistic particles and a wide range 
of matter densities (that can be found in diverse astrophysical and cosmological environments) 
the energy range of the SLv and SLe may even extend up to energies peculiar to the spectrum of 
gamma-rays (see also [1,3]). 

For the rate of the SLv in the case of ultra-relativistic neutrinos (p ^> to) we obtained [1,2] 

^SLu = 4^ 2 a 2 m 2 ,p, m v /p < a < p/m„, (69) 

where the matter density parameter a is given by (|15[) . in case of negative a the SLv can be emitted 
by antineutrino. The main properties of the SLv investigated in [1,2,39] can be summarized as 
follows [42]: 1) a neutrino with nonzero mass and magnetic moment when moving in dense matter 
can emit spin light; 2) in general, SLv in matter is due to the dependence of the neutrino dispersion 
relation in matter on the neutrino helicity; 3) the SLv radiation rate and power depend on the 
neutrino magnetic moment and energy, and also on the matter density; 4) the matter density 
parameter a, that depends on the type of neutrino and matter composition, can be negative; 
therefore the types of initial and final neutrino (and antineutrino) states, conversion between 
which can effectively produce the SLv radiation, are determined by the matter composition; 5) 
the SLv in matter leads to the neutrino-spin polarization effect; depending on the type of the initial 
neutrino (or antineutrino) and matter composition the negative-helicity relativistic neutrino (the 
left-handed neutrino vl) is converted to the positive- helicity neutrino (the right-handed neutrino 
vr) or vice versa; 6) the obtained expressions for the SLv radiation rate and power exhibit 
non-trivial dependence on the density of matter and on the initial neutrino energy; the SLv 
radiation rate and power are proportional to the neutrino magnetic moment squared which is, 
in general, a small value and also on the neutrino energy, that is why the radiation discussed 
can be effectively produced only in the case of ultra-relativistic neutrinos; 7) for a wide range 
of matter densities the radiation is beamed along the neutrino momentum, however the actual 
shape of the radiation spatial distribution may vary from projector-like to cap-like, depending 
on the neutrino momentum-to-mass ratio and the matter density; 8) in a wide range of matter 
densities the SLv radiation is characterized by total circular polarization; 9) the emitted photon 
energy is also essentially dependent on the neutrino energy and matter density; in particular, 
in the most interesting for possible astrophysical and cosmology applications case of ultra-high 
energy neutrinos, the average energy of the SLv photons is one third of the neutrino momentum. 
Considering the listed above properties of the SLv in matter, we argue that this radiation can 
be produced by high-energy neutrinos propagating in different astrophysical and cosmological 
environments. 

A remark on the possibility for Majorana neutrino to emit the spin light in matter should 
be made. Obviously, due to the absence of the magnetic moment, such radiation is not expected 
in this case. However, considering the transition between two neutrinos of different flavour, it 
is possible to produce an analogous effect via the transition magnetic moment, that Majorana 
neutrinos can posses. 

Performing the detailed study of the SLe in neutron matter [43] we have found for the total 
rate 

r SLe = e 2 m 2 /(2p) [In (4a„p/m e ) - 3/2] , m e /p -C a„ < p/m e , (70) 

where it is supposed that In °" p ^> 1. It was also found that the relativistic electrons can loose 
nearly the whole of its initial energy due to the SLe mechanism. 

It should be noted that discussion on possible impact of the background plasma on the SLv 
radiation mechanism have been started in [2] . Then effects of plasma for the SLv and SLe were 
considered by another authors in [44]. These authors, after we explained in [42] that their initial 
conclusion [45] that in presence of matter the process "r/j, — ► vr + 7* is kinematically forbidden" 
was wrong, obtained the SLv rate with account for the photon dispersion in plasma. In the case 
of ultra- high energy neutrino (i.e., in the only case when the time scale of the process can be 
much less than the age of the Universe) the SLv rate of [44] exactly reproduces our result (|69p 
obtained in [1,2]. The final result for the SLe total rate in the second paper of [44] in the leading 
logarithmic term confirms our result (|70[) obtained in [43] . 
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6. Conclusion 

We have considered a framework for treating different interactions of neutrinos and electrons 
in the presence of matter. The method developed is based on use of exact solutions of modified 
Dirac equations that include correspondent matter potentials. It has been demonstrated how 
this method work in consideration of different quantum processes that can proceed in presence of 
matter. 

Finally, let us consider the established in Sections 2.4 and 2.5 analogy between particles 
dynamics in the presence of electromagnetic fields and dynamics in matter. The developed 
semiclassical approach to description of the matter effect, driven by (electro)weak forces, is valid as 
long as interactions of particles with the background is coherent. This condition is satisfied when 
a macroscopic amount of the background particles are confined within the scale of a neutrino or 
electron de Broglie wave length. For the relativistic neutrinos or electrons the following condition 
should be satisfied the following condition must be satisfied - i 3> 1, where n is the number 

density of matter, 7/ = — L and (I = v or e). In case of varying density of the background 
matter, there is an additional condition for applicability of the approach developed (see, for 
instance, [19,27,46]). The variation scale of matter density should be much larger than the 
de Broglie wavelength, \-^-\ <C 1. 

We can further develop the established in Section 2.5 analogy between a neutrino motion in a 
rotating matter and an electron motion in a magnetic field. It is possible to explain the neutrino 
quasiclassical circular orbits as a result of action of the attractive central force, 

F&> = q% ] P x B m , B m = V x A m , A m = nv, (71) 

where the effective neutrino "charge" in matter (composed of neutrons in the discussed case) is 
qm = —G, whereas B m and A m play the roles of effective "magnetic" field and the correspondent 
"vector potential" . Like the magnetic part of the Lorentz force, F„ is orthogonal to the neutrino 
speed [3. 

It is possible to generalize the discussed above description of the matter effect on neutrinos 
for the case when the matter density n is not constant. For the most general case the "matter 
induced Lorentz force" is given by 

F&> = q^E m + gM/3 x B m , (72) 

where the effective "electric" and "magnetic" fields are respectively, 

^ _ dn dv , . 

E m = -Vn-v--n-, (73) 

and 

B m = nV x v - v x Vn. (74) 

Using ((4]) and (|5|) (see also (fT5|) ) these expressions can be generalized for a background composed of 
different matter species. The force acting on a neutrino, produced by the first term of the effective 
"electric" field in the neutron matter, was considered in [19]. Note that the same quasiclassical 
treatment of a neutrino motion in the electron plasma was considered in [46] . 

To conclude, we ague that it is also possible to introduce the "matter induced Lorentz force" 
for an electron moving in background matter. The weak forces acting on a neutrino and an electron 
in matter are identical. Therefore, similar to the case of neutrino, we can write for the force acting 
on an electron F„ in background matter 

F^?l e) E m + ? ^xB mi (75) 

where appropriate magnitude for the effective electron "charge" in matter qm should be used. 
As it follows from (|73[) and (fTS")) , an accelerating force acts on an electron when it moves in 
background matter with nonvanishing gradient of density. Using this observation, we should like 
to discuss a new mechanism of electromagnetic radiation by an electron moving in the neutrino 
background (m = v) with non-zero gradient of its density. This situation can be realized in 
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different astrophysical and cosmology settings. For instance, this phenomenon can exist when an 
electron propagates in the radial direction from a compact star object inside a dense environment 
composed predominantly of neutrinos that also move in radial direction after they were emitted 
from a central part of the star. In this case that total power of the radiation (in the quasiclassical 
limit) is given by 






,2 („a\2 



(a/3) 2 



.(l-/? 2 ) 2 ' (1-/3 2 ) 3 -' (?6) 

where (3 is the electron speed and a is the electron acceleration induced by the gradient of the 
neutrino background density. We expect that the proposed mechanism of the electromagnetic 
radiation can be important in other astrophysics settings like one that can be realized in neutron 
stars, gamma-ray bursts and black holes. 
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